www.gusucode.com > 小波分析理论与Matlab 7实现 > 小波分析理论与Matlab 7实现/code/Ch17/例程17-2.m
%装载语音信号 N=1024; s=wavread('m?l.wav',N); figure(1); plot(1:N,s,'LineWidth',2); xlabel('时间 n'); ylabel('幅值 A'); %选用小波函数 wavelet='db5'; %分解级数 level=5; %压缩参数 alpha=1.5; %阈值类型 sorh='h'; %小波分解 [c,l]=wavedec(s,level,wavelet); %计算阈值 [thr,nkeep]=wdcbm(c,l,alpha); %使用硬阈值压缩信号 [xc,cxc,lxc,perf0,perfl2]=wdencmp('lvd',c,l,wavelet,level,thr,sorh); %对系数编码 cC=encode(cxc); %传输的系数 tC=cC; %传输的系数向量长度 tL=length(tC); %零系数的百分比 pzeros=perf0 %保留的能量大小 pnormen=perfl2 %编码的压缩比 compratio=length(s)/length(tC)